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An explanation is proposed for the fact that Lepage-Mackenzie tadpole improvement does not work well for 
staggered fermions. The idea appears to work for all renormalization constants which appear in the staggered 
fermion self-energy. Wilson fermions are also discussed. 
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One-loop renormalization constants for a num- 
ber of staggered-fermion operators have large fi- 
nite parts (see e.g. ref. Q]), which appear not 
to be explained by gluon-tadpole contributions. 
This is in contrast to the situation with Wilson 
fermions, where tadpole improvement Q appears 
to work well. 

Staggered and Wilson fermions differ in the 
way they deal with the doubling problem: in 
the Wilson case, they get a mass of the order 
of the cutoff and decouple, while they are present 
in the staggered case, yielding a number of con- 
tinuum flavors which is a multiple of four. The 
staggered doublers do contribute to loop dia- 
grams, and here I suggest that this may explain 
why tadpole improvement does not work well for 
staggered fermions, by considering the one-loop 
fermion self-energy. 

Let us first consider the Wilson case. The one- 
loop self-energy is || 
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with eg i 2 given in Table |l| (in Feynman gauge; 
1st column). The 2nd column gives the tadpole- 
improved values; the 3rd the contribution to those 
of the 2nd from the integration region — tt/2 < 
ifj, < tt/2 (£ is the loop momentum). 

The table shows that gluon-tadpole improve- 
ment works well (I used the mean link in Feynman 
gauge). If we split up the Brillouin zone (BZ) for I 
as I = TTA+i, A = 1, . . . , 16, with I G (-tt/2, tt/2] 
and tta e {(0,0,0,0), (n, 0,0, 0), . . .}, we see that 



most of the tadpole-improved values comes from 
the region with tta = 0. The doublers do not con- 
tribute much, since they are suppressed by the 
Wilson mass term. 

Now, let us consider staggered fermions. The 
most general four-flavor mass matrix is 

M = m s +ml^ + -ml v {-i^ v ) 

+ mf £ 5 , (2) 

where the £ M are a set of gamma matrices in flavor 
space (m^ lv is antisymmetric). The labels S, V, 
etc. (for scalar, vector, etc.) denote irreducible 
representations of the staggered fermion symme- 
try group, and correspond to zero-, one-, etc. link 
operators in the lattice action. 

The staggered one-loop self-energy is 
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°t m (-i^ v ) + a A m 5 J^ 5 



-crp mf £ 5 



with U diagonalizing M: M d = UMW. The 
values of r and the ct's are given in Table || (ay 
is a similar constant for a one-link mass term for 
"reduced" staggered fermions Q ) . 

For the staggered case we see that gluon- 
tadpole improvement does not work for as, ay 
and ay (it does work for ar, a a and ap), 
and that the tta — region does not give the 
main contribution, giving a value typically much 
smaller than the tadpole-improved value. 
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no improvement 


tadpole improved 


reduced BZ 




0.326 


0.0158 


0.0169 




-0.0878 


-0.0103 


-0.00839* 


02 


-0.0120 


-0.0120 


-0.0184 


Table 1 








Wilson contact terms. * Depends on routing of external momentum, which can make a r 


^25% difference. 




no improvement 


tadpole improved 


reduced BZ 


T 


-0.0446 


0.0329 


0.0149* 




0.197 


0.197 


0.0276 


a v 


0.00385 


0.0813 


0.0244 


(T T 


-0.117 


0.0380 


0.0111 


0,1 


-0.209 


0.0233 


0.00516 


CTp 


-0.294 


0.0161 


0.000294 


OV 


0.0813 


0.159 


0.0332 


Table 2 








Staggered contact terms. * See Table ||. 







The interpretation I would like to put forward 
is as follows. Consider the standard one-loop self- 
energy diagram: 




For ir a 7^ 0, the gluon propagator is very sup- 
pressed (of order a 2 ), and the gluon line in the di- 
agram is effectively reduced to a four-fermion cou- 
pling between the staggered flavors. The fermion 
propagator, however, has poles for tta ^ 0, and 
the diagram contributes "doubler-tadpoles" to 
the finite part of the one-loop self-energy. If we 
contract the gluon line to a point, the integrand 
near the doubler poles (tta ^ 0) goes like l/£ 2 , 
and produces contributions very much like the 
usual gluon tadpoles: the integral is quadratically 
divergent (1/a 2 ), with an extra factor ~ a 2 com- 
ing from the contracted gluon propagator. 

This idea can be tested by replacing the gluon 



propagator 
gauge) by 
1 



(4£ sm 2 (V2) (in Feynman 
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±[p-(e 1 ) P -(i 2 )P+(i 3 )P+(e 4 ) + 
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16 
where 



[P-(h)P-(l2)P-(l 3 )P + (l4) + 
P-(h)P-(l 2 )P-(h)P-(h) , 



~(l±cos^) . 



(5) 



The factors p-^ 1 )P+(£ 2 )P+(4)P+(£ 4 ), etc. are 
"smooth projectors" onto the regions tta +^ with 
tta 7^ of the Brillouin zone. The smoothness 
makes it possible to interpret this "gluon" ex- 
change as a local four-fermion operator. The frac- 
tions 1/4, 1/8, etc. are the values of the real gluon 
propagator at I — tta ^ 0. Note that the region 
around tta — is "projected" onto 0. 

With this replacement, i.e. with this four- 
fermion interaction, one obtains the results of Ta- 
ble §. We note that 

• the four-fermion constants reproduce the 
tadpole-improved contact terms quite well, 
especially the larger ones; 

• gluon-tadpole improvement is also needed; 

• for Wilson fermions, the "four-fermion val- 
ues" for do, a i and a 2 are —0.00222, 
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tadpole improved 


four-fermion 


fraction 


T 


0.0329 


0.0185 


0.56 


CS 


0.197 


0.179 


0.91 


oy 


0.0813 


0.0668 


0.82 


(T t 


0.0380 


0.0309 


0.81 


OA 


0.0233 


0.0237 


1.0 


Up 


0.0161 


0.0227 


1.4 


a v 


0.159 


0.134 


0.84 



Table 3 



Staggered contact terms - comparison with four-fermion values. 



-0.00111 and 0.00949, respectively, to be 
compared with the gluon-tadpole improved 
values given in Table 0. 

The idea presented here can be checked on 
the many other staggered-fermion renormaliza- 
tion constants that have been calculated to one 
loop in perturbation theory. Also, since the ap- 
proach is gauge-dependent, it should be checked 
for other gauges, such as Landau gauge. 

Assuming that the idea is correct, one can ask 
how better estimates of staggered-fermion renor- 
malization constants can be obtained. This could 
be done by taking the staggered-fermion theory 
with only the four-fermion interactions (no glu- 
ons), and computing the renormalization con- 
stants in this theory numerically. (These con- 
stants are finite, because the four-fermion inter- 
actions are irrelevant operators, proportional to 

<r.i 

To one-loop, one then multiplies these by the 
perturbatively calculated constants of the full 
theory (with gluons), with the four-fermion part 
taken out, and the appropriate power of the mean 
link for gluon-tadpole improvement. For exam- 
ple, for the wave-function renormalization, we 
would get 



Z 2 = u^Zf 



3 y 



-JLlog« M - 



,(6) 



where uq is the mean link, Zf is the wave- 
function renormalization of the four-fermion the- 
ory, and At = 0.0144 is the difference between 
the first two numbers of Table ||. 

This procedure resembles gluon-tadpole im- 
provement, in that it partially resums the pertur- 
bative expansion, and it is equally heuristic. A 



disadvantage is that a numerical computation is 
needed in the four-fermion theory for each oper- 
ator. A complete nonperturbative determination 
m the full theory [| may therefore be preferable 
not only in principle, but also in practice. 

A different approach would be to consider im- 
proved actions for lattice QCD with staggered 
fermions (see e.g. ref. Q). For improved actions, 
one expects the couplings of high-momentum glu- 
ons and fermions to be smaller than in the unim- 
proved case, which would presumably lead to 
smaller finite parts of the renormalization con- 
stants. 

I would like to thank Jim Hetrick, Mike Ogilvie, 
Steve Sharpe and Doug Toussaint for discussions. 
This work is supported in part by the US Dept. 
of Energy under the Outstanding Junior Investi- 
gator program. 

REFERENCES 

1. A. Patel and S. Sharpe, Nucl. Phys. B395 
(1995) 701. 

2. P. Lepage and P. Mackenzie, Phys. Rev. D48 
(1993) 2250. 

3. L. Karsten and J. Smit, Nucl. Phys. B183 
(1981) 103. 

4. M. Golterman and J. Smit, Nucl. Phys. B245 
(1984) 61. 

5. N. Ishizuka, for JLQCD, these proceedings. 

6. K. Orginos and D. Toussaint, these proceed- 
ings. 



